NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



MASA 

Technical Memorandum 82074 

(NASA-TM-U2074) SNOM WATIiK baUIVALEM NU1-204!#0 

UETEmilKA'lION BY MICKOWAVb KABXOMfitKY (HAS A) 

22 p UC AJ2/Mi A01 CSCE 08L 

Uacld£i 

aJ/4J 1S349 

SNOWWATER EQUIVALENT 
DETERMINATION BY MICROWAVE 
RADIOMETRY 


A. T. C. Chang, J. L. Foster 

D. K. Hall, A. Rango and B. K. Hartline 


Hydrological Sciences Branch 


JANUARY 1981 


National Aeronautics and 
Space Administration 

Goddard Space Flight Center 

Greenbelt, Maryland 20771 



1 

j 



vSNOW WA I'I'R l-OUIVAl FNT Dl'TI RMINA I'lON 
HV MK’KOWAVl- KADtOMF I'RY 


A. T. C. ('hang. J. I . I'oslor. I). K. Hall, 
A. Kangu aiul B. K. Hart line 
Hydrologioal Sciences Brunch 


CODDARO SPA('F |■U(:HT Cl'N'I'HR 


(iRHFNBlvl r, MARYLAND 20771 



SNOW WATER EQUIVALENi DETERMINATION 
BY MIC ROWAVE RADIOMETRY 

ABSTRACT 

Snow Water liquivalent (SWII) is one of the most important parameters for accurate snowmelt 
rimolT prediction. C’onventionally, SWE is monitored using observations made at widely scattered 
point.s in or around specific watersheds. Remote sensors, which provide data witli better spatial 
and temporal coveraii; , can be used to improve the SWE estimates, Microwave radiation, which 
can penet “ate through a snowpack, may be used Jo infer the SWE, 

C'aleulations made from a microscopic scattering model are used to simulate the effect of 
varying SWE on the microwave brightness temperature, Data obtained from truck mounted, air- 
borne and space-borne systems from various test sites have been studied. The simulated SWE com- 
IKire.s favorable with the measured SWE. In addition, whether the underlying soil is frozen or 
thawed can be discriminated succe.ssfully on the basis of the polarization ofthc microwave radiation. 
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SNOWWATER I QUIVALHNr DHTI'RMINATION 
BY M /( ROWAVl: RADIOMKTRY 


I. INlRODUnUJN 

Snow iiocumulalion aiu! Uepletion is a highly critical parameter in the western United Slates. 
Meltwater rroiii imnsnt.iin snowpacks provides tnuch needed water for hydropower generation and 
irrigation. As much as 70 percent of the water supply for the wesiern Utiited Stales is directly 
derived from the spring runoff of the snowpack, Knowledge of the available water volume in arl- 
vance of the runnoff period will permit better managenient of the expected runoff in the spring. 

I he sjrring meltwater also provides re<;harge to the upper soil mantle for sprirrg wheat and other 
grains in agricultural regions, Information on the water equivalent of these non-mour'ilainous snow- 
inrcks eoidd indicate areas of potential moisture deficiencies or surplus. 

In order to monitor the snow water equivalent accurately for these hydrological applications, 
ineasurcinenls of snow covered area, snow depth, density and liquid water content are essential. 

In iuldilion to these snow properties, Ihe condition of the underlying soil is important for ostimal» 
ing (he aimnml of snowmelt water that will be retained in the soil mantle, (’onvontionally, these 
observations are eolleeled manually by snow surveyors using i?ki,s or snowmobiles, and aulomatieally 
from a lew unattended, i.solated stations which are instrumented with pro.ssure pillows and other 
in-silu sensors. I’hese data collection melliods are time consuming and measnremonls are taken 


at only a liniiUnl number of points along selected snow survey eonrses, ('onsequenlly, the .snow 


water equivalent estimate, s may differ significantly from the actual water equivalent of a snowpack 


because of the sparsity of the observations made in space awl time. The use of remote .sensing 
U'chntcjuos may offer a way to augment or complement the conventional observations by providing 


high spatial density and repetitive ob.servalions over entire wntersheds. 

Sensors on board present operational satellites scan the earth in the visible, near infrared and 


tliermal infrared portions of the electrornagnetie spectrum. Snow covered area estimates from 



plitnarjiy vi.sil'V images (or sovcial (est walersheds (lave been fmmil to be well correlated with 
the aeliiiti snowmelt nmot't'tRango ami Peterson, 1980). However none ot'these short wavelength 
sensors can |<rovkle (he information on (he depth and water equivalent of the snowpack that could 
fuiiliei imi'rove water resources management. Microwaves liave tlie capability of penetrating the 
snowpaek, thus that portion of the electromagnetic spectrum has considerable |H>tential for moni- 
toring snowpaels water equivalent and wetness ( Rango, et al., 1979), In addition, microwave 
radiation i)f about cm wavelength or longer may penetrate through a typical snowpaek (-^ I m 
depth) to provitic inlbrmation on ’oil cond it ion.s (frozen or thawed) beneath the snowpaek. These 
eoiulit ions greatly influence the predietkm of snowmelt runoff yields. 

in or del to qit.mtitalively determine the v ater equivalent of the snowpaek by microwave 
railiometiN . it is necessary to understand tlie betiavior of microwave radiation within the snowpaek. 

} his (Kiperdesu ilH's a model of die nuctrrwave emission characteristics of snowpacks for various 
snow comiitions. Microwave radiometric data measured by truck-mounted, airborne and space- 
boine sensor systems have been compared with the calculated results of a microscopic scattering 
model (('lung, et al.. 19?o), fhe snow depth and waiter e(|uiva lent estimates derived from the model 
calculations compare favorably W'ith the limited ground truth information. 

II. nilOlU 

The intensity of microwave radiation emitted from a snowpaek depends on the physical 
temperature, grain si/e, density and the underlying surface conditions of the snowpaek. Hy know- 
ing: these parameters, the ladiation emerging from a .snowpaek can be derived by solving the 
radiative tran.sfer equation. The radiative transfer equation for an a.\ially symmelric inhomogeneous 
medium can be written in the form of an integro-differential equation 

h — ~ o( ,v) |(x. /j) + n( x) >1 1 I - oH.x)) B(.,x) + 'bto( x) | P(x, p') Kx. (i') d/a' , ( 1 ) 

-I 

w here the radiation inten.sUy l( x, /a) is at depth x traveling in the direction making an angle who.se 
cosine is f.t with the normal toward the direction of increasing x (Figure I). The functions (i(.x). 



lt(x) iiiui l*( \ uie pivscrihcil functions of their arguments. They are referred to as the 
extinction per unit length, the single scattering albedo, the senirce and the phase* functions, respec- 
tively . file snow grains scatter the electromagnetic radiation incrrhercntly aiul are assiiiiKnl to be 
spherical in shape ami raiuloinly spaced within the snowpack. Further discusshms of the micro- 
scopic model can be found in t'hung et al., (l‘)7(>). Fquation (1) is solved numerically by 
the invariant imbedding technique ((’hang and t’houdhury. I‘)78). This technique is based on the 
principle that the radiation emerging from a semi-infinite, plane pavalle! medium is invariant with 
respect to addith>n(or subtraction) cf layers of arbitrary thickness to (or from) the medium. 

I xpressed numerically . this technique solves the radiative transfer equation via recurrence relations, 

Hy ihis method, the emerging brightness temperatures are calculated for snowpacks with 
difleient parameters. I he siurvv depth varies from 10 cm to 250 cm while the water equivalent 
varies liom 5 cm to 100 cm. fhe mean snowpack grain si/c varies from .25 mm to .5 mm in railius. 

1 w o vliirerent ty pes of undeiiy ing surface conditions are considered, They are ( I) fro/en srri! and 
(,’) imiro/en soil w ith 20'i moisture content. Figures 2 and 3 show the calculated brightne.ss tem- 
perature as a luKction id' siu»w water ei)uivalent (SWF) for 37 (ill/ radiation with a look angle of 
50 degiees. In this paper we avhlre.ss rmly the ilry snow cinulition; the wet .snow condition will be 
aililresseil in a .sub.sei|uen( paper, 

Ua.si l i>n these calculaturns we may conchulc that the brightne.ss temperature at M (>11/ tends 
ti> decrease rapidly with inerea.sing SWF. AlSi», the calculated rate of change of brightne.ss depends 
srrongly on the grain si/e, yet it is almo.st independent of the siiow den.sity, This character i.stic may 
hi* u.sed as the basis to ileterntine the SWF. llow'cver, the emerging 37 (ill/ brightne.ss from a sm>w- 
pack depends i>n the underlying .soil conditions. In order to rest>lve this ambiguity, mea.surements 
frr>m another frequency preferably longer wavelength are required, five mea.surements at 10,7 (ill/ 
piarvide sufficient information needed for (his purposr'f. 
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Water u* a polar nu'lecuk and it has a very large dielectric conslaitt at itiicrowave frequencies. 
I'his results in a large rellectivity difference from a frozen and an unfrozen soil st;rface. For inci- 
dence angles not equal to zen>, the polarization factor 


is commonly used to reneci the different surface conditions. The Cit/Ctdaled polarization factor 
for 10 7 (lllz ranges Irom 0 OH to 0.1 1 S for unfrozen sod with 'O'i sv^>il moisture and from O.C '2 to 
0 051 lor liozen soil, riuishy a simple dt.sci itninant analysis, these two categories can he separated 
h> using the polarization factor: 


T > 0.07 

t atiweii soil 

i.D 

T w 0 07 

tio'cn soil 

(4) 


111. COMl’ARISON.S 

l licvc haw been several dilTeren! expenments designed lo ineasun the inieiowave Inightness 
teinpei ituie ,is .i lunetiou olMuns de|i(h or vsatet eiiuivalent. The first truck niounteil radioiiKier 
experiment was eoiuliieteil on a site m the northern part ot the Sierra Nevada Range near I tuekee. 
C'A by I dgertoii.et al,, l‘>73, I'hey reported a smooth and pionounced reduction ot 57 (.ill/ bright- 
ness temperatuie during (he snow accumulation period, rtiese results iiulicate that microwave 
brightness temperatures may I'c used to measure the water equivalent ofa snowpaek. Recently, 
several groups have been involved in microwiwe .snowpaek characterization experiments using 
ground-luised .sy.stems (Stiles aiul Ulaby, 1^80 llofer and Matzler, l‘>80, and ('nang, et al.. ld?d). 

I'he results reported by various investigators are consistent. Due to different snow conditions 
in different geograplucal locations, however, some di.ssimilarities have ah . j been reported. 
I'ypical brightness temper, itures at lillz and the related snow information from ilitfeient experi- 

ments ate tabulated in fable I Tor coimxirison purposes, the predicted SWT. is obtained by 
converting (lie measured brightness using the curves in Tigiires 2 and .T THie to the lack 
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of ailcqiute fucusureinenJ!* of t(>e mean snow crystal radius in each experiment , only one mean 
radius value (0.3S mm) was assumed for this initial simu’ition of the SWP. in all test sites. 

from fable I H istpiite encouraging that the predicted SWF. values are compaitble with the 
measuretl results for all the snowpacks. except in the deep SWF in Switzerland. Generally predic- 
tions are successful for snowpacks shallower than one meter. At Frazer. CO the percent polariza- 
tion faett)! ranges from 0.02 to 0,04. These values corres|H>nd to the frozen soil condition as 
obsem'd at this particular valley site. 

Snow crystal size usually increases as measurements are made deeper In the snowpack ((Jow, 
|d(e>), rhe radius for new dry powder snow on top of a snowpack normally ranges from 0,05 to 
0,1 mm, Whereas depth hoar, located at t ie bottom of a snowpack, may have a radius from I nun 
to several mm. Since at preseiU, the mean crystal size is assumed to be a fixed value (0.35 mm) In 
this study, It is possible to adjust the mean snow crystal radius as a parameter to obtain a better fit 
ol the present data. l‘or example. If we choose the radius to be 0.3 n.m, then tlie brightness tem- 
perature observeil at Fraser iiuUeatesa SWF of 25 cm instead of 16 cm as obtained using 0.35 mm 
for mean rail Ins. riierefore, in order to obtain reliable estimates of SWF by this lechnii|ue, it is 
necessary to eareiuliy characterize the physical size of the i rystals within the .snowpack. 

Use of tlie 3? (illz brightnesis to predict the SWF is limited, because the 37 Ghz radiation can 
only penetrate (fte top 50 to 100 cm of tl.e snowpack (Stiles and Ulaby, 1980. and Holer and 
Mat/ler, 1979 ), I'hus in case of a snowpack several meters in depth, the mean crys(.d radiiis should 
only he calculated by usiitg the crystal size from the uppermost portion that the microwaves sample. 
This top layer probably has a much smaller mean radius than the entire snowpack, Recent measure- 
ments by Matzler et al. (1980) in Switzerland showed that the mean radius of the 1979 snowcover 
is alioui 0.25 mm. If this value is taken as the crystal radius, then the SWF predicted using Figure 
closely apinoximates to the measured SWF of 50 cm. 
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Dunng tlK* winters ol 1976, 1977 and 1980 several aircraft mtssiotis took place in the Colorado 
Rockies for snow pack studies Steamboat Springs Rabbit F.ars Paw and Waklen.CO were the major 
study sites, The radiomeieis on board the aircraft include a four channel MuliM'requency Micro- 
wave Radiometer (Ml'MR) and a Passive Microwave Imaging System (PMIS). C oordinated with 
these flights, exlei dve ground truth information was gathered. Snow density, water equivalent, 
depth, wetness, appii'-ximate gram ta/e, lay*‘r classification and soil condition were measuicd by the 
ground truth teamst Jones, 1976. |977, 1070, 1980a and 1980b). For the puiposc of this study 
the ilata were averageil o\ct the eniiie flight line for better comparistm with the large microwave 
radiometer liiotprmt obtained froir. airborne sensors. The ,J7 (Jll/ brightness is used to predict tlie 
SWF while (he 10. till/ biighlness is im‘d to determine the underlying soil conditions, Table 3 
shows the results of neasured SWb as compared wUh SWF derived from microwave brightness I he 
comparison between the piedietv'd and measured SWi* ,s good everywhere except ai Rabbit bars 
Pass. Again, the diserei anev probably is due the arbitrary choice of 0.35 mm for the assumed 
grain si/e used m the model. ,\t Rabbit bars Pass tlie snow was about 1.5 m m depth, I ven w ith 
such di'cp .snow , the percent polar;,/,i!ion lactor derived liom the 10.7 ril!/ measurement does 
ivfleci the underlying soil c»>iulitions when using 0.07 as the divider between the unfro/en and 
fro/en soil, 

Microwave measurements from spaceborne systems have been available since tlie launch of the 
Nimbus-5 Stitellite in December. 1972, Rango, et al. ( 1976) and f oster, ct al. ( |980i used Ihe 
Nimbus 5 and 6 I leclneally .Scanned ‘uieiowave Railiometer ( l-SMR) data In an attempt to deter- 
mine if a eurielation existed between snow deptii aiul brightness tcmpeiature in relatively homo- 
geneous regions of (he United States, C'anada and Russia. Statistically significant regression 
rclalk'iisliips and reasonably liigli coefficients of determination (R** "O-S) were obtained in these 
studies Tliese i\ suits demonstrate a potential for estimating snow depth by passive microwave 
data from spaceboine .sensors, 



In thic study we will report loinj pr«liniinRry resuHf derived from the Nimbus*? Scanning 
Multichannel Microwave Radiometer (SMMP.'L SMMR is a five frequency, due! polarkeu micro- 
wave radiometer which measures the upwethng microwave radiation at 6,6, 10,7, 18,0, 2 1 ,0, and 
37,0 (]ll/ while scanning 25** to either side of the spacecraft with a constant incidence angle of 
approximately 50“ with respect to the Earth’s surface, The spatial resolution varies from 25 km 
for the 37 CIHz to 150 km for the 6,6 GH^, l^ctailed descriptions of this instrument can be found 
in the Nlmbus-7 Uscr'sGuide (Gloersen and Hardis, 1978), 

The study areas chosen were (I) central Russia and (2) the high pbins of Canada (Foster et al,, 
1980), The vegetation, to|K)graphy, climate and latitude of these two areas are similar. The gen- 
erally Oat terrains of these areas, which is sometimes broken by hills is covered with various grasses. 
Botli of these two areas experience very cold winters with snow possibly covering the ground from 
December to Mi!i' 'h, 

Due to limited available snow course data. It is rather difficult to compile SWK data In a timely 
fashion for comparison with the satellite measurements. In this study SWE is calculated by multi- 
plying the .mow depth and the snow density. Snow depth information was obtained from meteoro- 
logical stations, whereas the snow density was assumed to be 0,3 g/cm^ for all the calculations 
(variations from this assumed density could account for some of the scatter observed in the data). 
Figure 4 shows the 37 GHz vertical polarization brightness temperature versus the snow depth for 
thet'anadian and Russian test sites as compared with the calculated results. The figure shows that 
most of the data points fall in the range of mean radius of 0.3 to 0.5 mm. These results provide us 
more confidence in the assumed mean radius (0.35 mm) used in this study. The time period of the 
SMMR data used is February 15 to 2 1 , 1979, During this time period the snow depth variation 
across the study areas is quite limited (from I to 30 cm). Although the foot-print of each bright- 
ness measurement represents a rather large area (25 km X 25 km), the effect of snow cover can still 
be observed. Figure 5 shows the scattering of SMMR 37 GHz brightness temperature ver.sus snow 
depth for the Russian test site. The linear regression technique gives a R^ of 0,75 for the Tg and 


7 



snow ilepil). The data display considerabb scalier, which is probably due lo Ihe iiilioiiiogeneily 
wiihhi each fuoiprinl and Ihe assuinplioiis used in this sludy. The Iheoretioiliy calculaled bright’ 
ness curve fils well with observations and can be utilized to infer the snow depth or SWE from 
inicrowK brightness lemperatuie measured by spaceborne sensors. No attempt has been made in 
this study to relate the polarization factor to the free/e/thaw soil condition determination, because 
there is no reliable “ground truth inforfliation" for comparison. 

IV. (ONCI USIONS 

I he microwave radiometer obscrvation.s can be utilized to infer the snow water equivalent 
under dry snow conditions. The results from a scattering model match well with the experimental 
results. I he capability of the M CJIIz radiometer has been demonstrated by using data collected by 
ground experiments, aircraft and s|>acecraft measurements. For deeper snowpacks, a wavelength 
longer tlian O K cm is required to infer Ihe snowpack information, The polarization factor derived 
from Ihe 10,7 (ill/ brightness temperature pixrvides an adequate index of the underlying soil 
corn! it ion. 

Due to the .strong dependence of Ihe enierging brightness temperature on Ihe mean crystal 
radius, it is necessary lo carefully eharaeteri/e (he crystal sizes within the snowpack. In addition, 
crystal size distribution should be monitored and documented in order to account for the different 
type of snow metamorphism found indifferejit snow sites. 

At present, the poor sensor resolution from .siitellite observations limits the use of the isaiellite 
to large homogeneous regions such as Ihe high plains. Even this coarsely derived information could 
be valuable for runoff prediction purposes for the lime period right before the rapid spring melt in 
these test sites. In thc,se areas, snow also provides the necessary insulation for Ihe underlying vege- 
tatioii such as winter wheat. Satellite derived information on snow depth could be the key to early 
detection of w inter kill. This information will greatly enhance the accuracy of overall crop yield 
prediction. As satellite spatial resoluUon improves with, future generations of microwave radiom- 
eters, applications should be found in other areas suc!i as intermountain valleys and large mountain 
plateaus. Then the capabilities of microwave radiometry will become more directly applicable to 
seasonal and short term runoff forecasting, 
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Table I. Compariion of Trusit Mounted Experimental Data 


Test Site 

T 

incidence 

Polarization 

Snow Water E;.{Uivalent 


Angle 

Measured 

Predicted 

Truckee, CA 

233K 

45* 

V 

10 cm 

9 cm 

(Kilgerton, e! al,, 1973) 

210K 

45* 

V 

20 cm 

16 cm 

Steamboat Spring, CO 

20SK 

57* 

H 

10 cm 

9 cm 

(Stiles and Ulaby, 1980) 

I88K 

57* 

H 

20 cm 

15 cm 

Fraser, CO 
(( hang.et al.. 1979) 

210K 

O 

O 

V 

21 cm 

16 cm 

Ddvus, Switzerland 
(lloferand Mat/.lcr, 1980) 

210K 

f 

V 

50 cm 

16 cm 


II 




Table 2. ( omparison of Aircraft Expcrimniial Data 


Test Site 

>n 

^11 

“10.7 

Polar- 

ization 

Factor 

Snow Water Equivalent 

Cl round 


(3"'H) 

Measured 

PrtHlicted 

Cover 

Steamboat Springs, (X) 
March 

105 

245 V 
lOOH 

0.10 

lo.J 1 7cm 

17 cm 

Wet 

WaKleii,( () 
Maid) l‘)7o 

235 

2 08 V 
258 ii 

0,02 

2.8 i 1,5 cm 

3 cm 

Fro /on 

Steiimbuat Springs, (X) 
Match l‘)77 

200 



10.4 s 3,7 cm 

1 1 cm 

Fro/cn 

Wakicn.t O 
March l')-7 

254 



0.12 cm 

0 cm 

Fru/en 

Sii'aiiilu'at Springs. (X) 
iH'cemtu'. I‘)7‘) 

240 



0.0 1 O.o cm 

2 cm 

F'ro/eii 

Kabl^il bars Pass.ro 
l•Vl1IuaIy l')S() 

234 

2 03 V 
245 M 

0 04 

44.3 t 13 cm 

4 cm 

i'l'o/cn 

Sleamlnr.il .Springs. ( () 
I'ehiiiaiN iMSt) 

I'O 



14.8 1 2 cm 

1 2 cm 

l''ro/cn 

Sic.imhoal Springs, CO 
March I'hst) 

lo2 

255 V 

210 11 

U.U) 

30,4 t 5.0 cm 

24 cm 

Wei 

W'alilen, ro 
Marcli !‘)SU 

104 



4,4 1 4,8 cm 

1 2 cm 

Wet 



RGUKli CAmONS 


I igiire I . Kiiiitutioii Inlcnsily of l(x, m)- 

linurc 2 , (a leu 111 toil .17 (111/. BrislUncss Temperature vs. Snow Water Equivalent Over Frozen 
(Jrouiul (Incidence Angle * 50°), 

l iguiv .1. C alculated .17 (111/ Bright ness Temperature vs, Snow Water Equivalent Over Unfrozen 
Soil (Incidence Angle * 50°), 

I'igure 4. (’ompaiisons of ('ulculated Brightness vs, Measured Brightness for Different Snow Depths 
(,17 (iil/ Vertical I’oluri/ation. 0 = 50°). 

Figure 5, NIMIUJS-7 SMMK .17 (III/, Vertically Polarized Microwave Brightness remp,'iature ( I'lj) 
vs. Snow Depth (Russia) R“' = 0,75. 
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